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Abstract. Let {Mi}f =1 be a non-trivial family of d X d complex matrices, in the 
sense that for any n £ N, there exists i\ ■ ■ ■ i n £ {1, . . . , £} n such that ■ ■ • Mj n 7^ 
0. Let P: (0, 00) — > K be the pressure function of {Mj}f =1 . We show that for each 
g > 0, there are at most d ergodic ^-equilibrium states of P, and each of them 
satisfies certain Gibbs property. 



In this paper, we study the thermodynamic formalism for matrix products. We will 
characterize the structure of equilibrium states of pressure functions, and also examine 
the Gibbs properties of such states. This work was first carried out in [TT] in the case 
that the involved matrices are non-negative and satisfy a kind of irreducibility. Some 
applications were given in the multifractal analysis of the top Lyapunov exponents 
of matrix products [HI EJ [8] (see also [TO]). In this paper, we will consider arbitrary 
complex matrices. 

Let (£, a) be the one-sided full shift over the alphabet {1, . . . ,£} (cf. [2]) and let 
{Mi} e i=1 be a family of d X d complex matrices. For q > 0, we define 



where S n is the collection of all words of length n over {1, . . . , £}, Mj = Mj 1 ■ ■ ■ Mj n 
for J = ji • ■ ■ j n , and || • || is the standard matrix norm. By sub-additivity, the above 
limit exists and takes values in the set K U {—00}. The function P is called the 
pressure function of {Mi} e i=l . It plays an important role in the multifractal analysis 
of Lyapunov exponents of matrices [TTJ [6j [8] . Moreover, it is closely related to the 
dimension theory of self-affine sets and measures jH [15]. 



1. Introduction and results 
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Denote the collection of all cr-invariant Borel probability measures on £ by 
Endow jMo-(£) with the weak-star topology. For \i e A^o-(S), we define 

(1.2) MM = lim - V ^i([J]) log HMjII, 

n->oo z — ' 

Jes n 

where [J] denotes the n-th cylinder {x = ( = J} in E. The term 

M*(/i) is called the Lyapunov exponent of {Mi} e i=1 with respect to fj,. It also takes 
values in the set 1U {— oo}. The following variational principle for P was proved in 
[3] in a more general sub-additive setting: 

(1.3) P{q) = swp{qMM + h(p) : li E M a {Z)}, 

where h(n) denotes the measure-theoretic entropy of li with respect to a (cf. [19J). 
We remark that (11.31) was proved earlier in [3 [15] when the matrices are non-negative 
or invertible, respectively. For given q > 0, let 

(1.4) X q = {li e M a {E) : P{q) = gM*(//) + %)}. 

Each element li in X q is called a q- equilibrium state of P. Since both M*(-) and /i(-) are 
upper semi-continuous on A^ CT (S), X g is a non-empty closed convex subset of •Mo-(E). 
In particular, X q contains ergodic elements (each extreme point of X q is an ergodic 
measure). 

Our main purpose is to characterize the structure of X q . This question was par- 
tially raised from [16]. A complete characterization is given in Theorem 11.71 In the 
following, we shall present the setting and results. Proofs of the results are postponed 
until §2 

Definition 1.1. Let F be 1R or C. A family of d x d matrices {Mj}f =1 with entries 
in F is said to be irreducible over ¥ d if there is no non-zero proper linear subspace V 
of ¥ d such that M,V C V for all i £ {1, . . . , £}. 

The above definition is adopted from (TJ p. 48]. If {Mi} £ i=1 is irreducible over 
¥ d , then there exist D > and k G N such that for any words /, J G X* = 
U^Li{1j • • • > Q") there exists a word K in |J n=1 {l, . . . , £} n such that 

(1-5) WikjW >D||Mj||||Mj||. 

For a proof, see [B, Proposition 2.8]. This property is crucial in the proof of the 
following proposition. 
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Proposition 1.2. Let F be M or C, and {Mj}f =1 a family of d x d matrices with 
entries in F. If {Mj}f =1 irreducible over¥ d , then for each q > 0, P has a unique 
q-equilibrium state \x q . Furthermore, \x q has the following Gibbs property: 

(1.6) L7- 1 exp(-nP(g))||M J r <//,([ J]) < Cexp(-nP(g))||Mj||' 

/or all n £ N and J G £ n . Moreover, P is differentiable over (0, oo) and P'(q) = 
M*(n q ) forq>0. 

Remark 1.3. Proposition 11.21 is an analogue of Bowen's theory about the equilibrium 
state of Holder continuous additive potentials (cf. |2J). See jTH [19] for backgrounds 
and more details about the classical thermodynamic formalism of additive potentials. 
Proposition 11.21 was first proved in [TT] for non-negative matrices under a different 
irreducibility assumption (that is, there exists r G N so that X^: =1 (Mi + ' ' ' + 
is a strictly positive matrix). An extension was recently given in |9] Theorem 5.5] to 
certain sub-additive potentials. 

Let us next consider the non-irreducibility case. Denote the n x m zero matrix by 

Onxm' 

Proposition 1.4. Let F be K or C, and {Mi} e i=1 a family of d x d matrices with 
entries in F. Then there exist an invertible d x d matrix T, t G {l,...,d}, and 
positive integers d±, . . . , d t with d = d\ + • • • + d t such that for every i G {1, . . . , £} the 
product T~ X M{T is a partitioned matrix of the form 

(1.7) T- X M{T = (A[ j ' k)S 



l<j,k<t ' 



where A^ ,k \ j, k G {1, . . . , t}, satisfy the following two properties: 

(i) A^f'^ is a dj x dp. matrix and A^' k ^ = Od jX d fc when j > k. 

(ii) For any j G {1, . . . , t], either the family {A^' 3 ^}f =1 is irreducible over F d ^, or 
A^ = djXdj for allie{l,. ..,£}. 

Considering the partition (11.71) in the above proposition, we set 

A = A({Mi}f =1 ) = {j G {1, . . . ,t} : {4 0j) }i=i is irreducible over ¥ d ^}. 

Remark 1.5. It is possible that A = 0. For instance, this is the case for {Mj}? =1 , 
where 

m 1 = i' 01 V m 2 =( 02 

1 / ' loo 
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Anyhow, it holds that A = if and only if there is k G N such that M i± ■ ■ ■ M in = dxd 
for all n > k and if--i n G {1, . . . ,£} n . Observe first that T~ l ■ ■ ■ M in T is a 
partitioned matrix of the form (-B^' fe ^)i<j,fc<t, where 

^ g , l<yi,..,y n -i<t 

i<?;i<?y2<"-<i/n-i<fe 



is a <ij x <ifc matrix. According to (ii) of Proposition 11.41 A = implies A^' = 
for alH G {1, . . . , £} and j G {1, . . . , t}. Hence M h ■ ■ ■ M in = dxd for all n > t by 
(II. 8p . To see the converse, assume contrarily that {A^'^Y i=l is irreducible over F dj ' 
for some j G {!,...,£}. It follows now from (11.51) that for every n6N there exists a 



dj x dj 



word n • • • i„ such that Af^' ■ ■ ■ Af^ ^ ®d 3 xd 3 and, consequently, M h ■ ■ ■ M in ^ d 



xd- 



Definition 1.6. A family of d x d complex matrices {Mi} e - =l is called non-trivial 
if A 7^ 0, or equivalently, for each n G N, there exists / G {l,...,£} n such that 
Mj ^ dxd . 

In the following, we always assume that {Mi} e i=l is non-trivial. If j G A, then 
the pressure function of }f=i is denoted by Pj and the Lyapunov exponent of 

{>V }f=i with respect to /i is denoted by A*(fi). The following is the main result 
of our paper. 

Theorem 1.7. In the above general setting, it holds that 

(i) M*(//) = max{yll^(/i) : j G A} /or eac/i ergodic measure /j G A / l cr (S). 

(ii) P is a real-valued convex function on (0, oo), andP(q) = max{P, (g) : j G A} 
/or all q > 0. 

(iii) if q > and fij tQ , j G A, z's £/ie unique q- equilibrium state for Pj, then 

X q = conv{/i Ji9 : P, (g) = P(g)}, 
where conv^A) is the convex hull of A. 

Remark 1.8. The equality in (i) of Theorem 11.71 may fail for non-ergodic measures of 
jMo-(E). For instance, consider {Mi}f =1 , where M x = diag(l,2) and M 2 = diag(3,2). 
Let (here 5 X denotes the Dirac measure at x), and fi = pfi\ + (1 — 

p)/i 2 for some < p < 1. It is easy to check that 

M*(//0 = log 2, = 0, 4 2 Vi) = log2 
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and 

M,M = log 3, A^(^) = log 3, 4 2 H/i 2 ) = log 2. 
Since M*(-), A^\-), and ^i 2 ^-) are affine on J\A a (E), we have 

MM =plog2 + (1 -p)log3, A^(fi) = (1 -p)log3, A£ 2 >(//) = log 2, 

and thus, M+{p) > max{A^(ji) : i G {1,2}}. 

Remark 1.9. The pressure function for products of matrices has been studied in the 
literature under some stronger conditions. Let {Mj}f =1 be a family of real invertible 
matrices. Assume that {Mj}| =1 satisfies the strong irreducibility and contraction 
conditions (cf. p],[T3]). Guivarc'h and Le Page showed in [131 Theorem 8.8] that the 
pressure function P of {Mj}f =1 corresponds to the logarithm of the spectral radius 
of certain Ruelle transfer operator and moreover, P is real analytic on (0, oo), and it 
can be extended to an analytic function on {z G C : $t.z > 0}. This strengthens an 
early result of Le Page [17] . 



2. Proofs of the results 



This section is dedicated to the proof of Theorem 11.71 For the convenience of the 
reader we shall also present complete proofs for Propositions 11.21 and 11.41 

Proof of Proposition II . ,81 Let q > 0. Define a sequence of probability measures 
{vnq)n>\ on S so that 

rrm W|! 
V " AI]) = E, EE „ W 

for all I G S n . Let v q be a limit point of the sequence {v nq ) n >i in the weak topology. 
Furthermore, let fi q be a limit point of the sequence 

n-l 



n 



v q o a 3 



n 

J=° ' n>l 



in the weak topology. Using (11. 5R and a proof essentially identical to that of [TT1 
Theorem 3.2], we see that \i q G A4 a (E) is ergodic and has the Gibbs property (II. 6p . 
Thus 



+ > lim - V ^([J])log(L7- 1 exp(nP(g))^([J])) 

JGS n 

- hm - £^([J])log^([J]) = P(g). 



n— >oo IT, 

Jes„ 



6 



DE-JUN FENG AND ANTTI KAENMAKI 



Recalling fll.3j) . this implies fi q 6l g . 

Applying ( ll.6p and the ergodicity of and using an identical argument as in [21 
proof of Theorem 1.22] (or using [161 Theorem 3.6]), we see that fi q is the unique 
element in X q . According to this uniqueness, we have P'(q) = M*(/i 9 ), which follows 
from the Ruelle-type derivative formula of pressures obtained in [7J Theorem 1.2]: 

P\q-) = mf{M*(fi) : /i G XJ, P'(q+) = sup{M*(/i) : fi G XJ. 

We remark that although [71 Theorem 1.2] only deals with non- negative matrices, the 
proof given there works for arbitrary matrices. Alternatively, to show that P'(q) = 
M*(/x g ), we may apply ( II. 6p and the ergodicity of fj, q , and follow [HI proof of Theorem 
2.1] (see also [151 Theorem 4.4]). □ 



Proof of Proposition \l-4\ We prove the proposition by induction on d. Clearly the 
proposition is true when d — 1. Assuming there exists an integer p so that the 
proposition is true for all d < p, we show below that it remains true for d = p + 1. 
Let L(n,m) be the collection of all n x m matrices with entries in F. 

If {Mi}f =1 is irreducible over ¥ d , we simply take t = 1 and have nothing else to 
prove. We may thus assume that {Mj}f =1 is reducible, that is, there exists a non- 
zero proper linear space V of W d such that M^V C V. If we let v = dim\^, then 
1 < v and d — v < d — 1 = p. We choose an invertible linear map T\ : ¥ d — > ¥ d 
such that Ti(¥ v x {0}) = V. Then for each i e {!,...,£} there exist E t G L(v,v), 
Bi G L{y , d — v), Di £ L(d — v,d — v) so that 



d)xd -D-i 

Now by the induction assumption, there exist invertible matrices T 2 G L(v,v) and 
T 3 G L(<i — t>, <i — v) such that (T 2 " 1 £ , jT 2 )f =1 and (T 3 " 1 X'jT3)f =1 have the desired 
partitioned form for all z G {1, . . . , £}. It follows that 



T 4 = T! 

is an invertible d x d matrix and 



T2 Ovx(d-v) 
O(d-u)xi) X 3 



1 , , I T 1 EjT<2 Try BjT-i 

T7 l MiT A = 2 1 

J(d-v)xv J 3 -^i-fS 



has the desired partitioned form for all z G {1, . . . , £}. □ 
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Before proving Theorem 11.71 we shall first prove the following auxiliary result. 

Proposition 2.1. Let (X, J 7 , /i) be a probability space and T : X — >■ X an ergodic 
measure-preserving transformation. Let {/ n }^i be a sequence of non-negative Borel 
measurable functions on X such that sup a . gX f\(x) < oo and 

(2.1) f n+m (x) < f m (x)f n {T m x) 

for all m,n G N and x G X. If e > and a = lim n ^. tX) (l/n) J \ogf n d/J,, then the 
following claims hold: 

(i) If a ^ —oo, then for [i- almost every x G X, there exists a positive integer 
n (x) such that 

(2.2) \\og f n (T m x) -na\ < (n + m)e 

for all n > n (x) and m G N. 

(ii) If a = —oo, then for any N > and fi-almost every x G X , there exists a 
positive integer no(x) such that 

(2.3) log f n (T m x) < -Nn + [n + m)e 
for all n > no(x) and m G N. 

Proof. We only prove (i). The proof of (ii) is similar. 

Assume that a G M.. Let e > and take < 5 < e/4. By the Kingman's sub- 
additive ergodic theorem, for /z-almost every x G X, there exists uq(x) such that 

| log f n {x) — na\ < n5 

for all n > n (x), and 

| log/ m (x) - ma | < (no(x) + m)5 
for all m G N. Hence by (12. ip . we have for n > Uq{x) and m G N, 
log/„(T m x) > log/ n+m (a;) - log/ m (x) 

(2.4) > (n + m)(a — 8) - m(a + 8) — n (x)S 

> na — 2(n + m)8 > na — (n + m)e. 



To see the opposite inequality, take k large enough such that |/3 — a| < 8, where 

1 
k 



P = T I lc, g fk d[i. 
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Since {fn( x )}%Li is sub-multiplicative, by [31 Lemma 2.2], we have for any n > 2k, 

n—k 

(fn(x)) k < C 2 " 2 J] f h (T*x) 
3=0 

for all x G X, where C = max{l, sup xgS f\(x)}. It follows that for n > 2k and m G N, 

n—k+m 1 m— 1 

log/.(T m x)<2fclogC+ £ -log/fc^-^-log/fc^). 

8=0 1=0 

Applying the Birkhoff ergodic theorem to the function \ log and combining it with 
the above inequality, we see that for //-almost every x G X, there exists an integer 
^o(^) > 2k8~ 1 \ogC such that 

log/ n (T"V) < n5 + (n - k + m)((3 + 5) - m(fi - 5) + h (x)5 

< n(3 + 2(n + m)8 + fio(x)5 < not + 4(n + m)<5 

< na + (n + m)e 

for all r/ > ho(x) and m G N. This together with ( 12. 4 j) yields ( 12.2ft . □ 

As a direct corollary of Proposition 12.11 we have the following. 

Corollary 2.2. Under the assumptions of Proposition IK7[ for any e,N > and for 
/i- almost every x G X , there is C(x) > such that 

|/ n (T m x)| < C{x) exp(n maxja, — iV}) exp((r/ + m)e) 

for all n, m G N. 

Proof of Theorem \1.7\ We only need to prove part (i), since parts (ii) and (iii) follow 
immediately from (i), the variational principle ( II ,3p . and Proposition 11.21 

Fix an ergodic measure // G M. a {T). The direction M+([i) > max{A* (//) : j G A} 
follows from the fact that 

|Up)... A^ j) \\ < llT^Mi, ■■■M l T\\< lir-^lllTllllMi, ■■■M i II 

1 1 %\ l n 1 1 — II L l L n II — II 1 1 1 1 1 1 1 1 L l L n 1 1 

for any j G A and ii, . . . , i n G {1, . . . , £}. We only need to prove the other direction. 

By Furstenberg-Kesten's theorem [12] on random matrices, or Kingman's sub- 
additive ergodic theorem (see e.g. |19j), we have for //-almost every x = (xj)^i G S, 

(2.5) lim - log || M Xl ... Xn || =M*(//). 
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For any i G {1, . . . , £}, define a sequence {fn }%Li of non-negative functions on E by 
setting 

fSHx) = \\AM...A<ip\\ 

for all x = (xi)^! G S. Let e, A > 0. Apply Corollary O for {fn j) }^ =1 to obtain 
that, for /i-almost every x = (x i )°l 1 G £, there exists C(x) > 1 such that 

, N U { ^l Xm+2 ... Xm+n \\ < C(x) exp(nmax{i4?)(Ai), -iV}) exp((n + m)e) 

(2.6) 

< C(x) exp(nmax{W, —A}) exp((n + m)e), 

for all j G {1, . . . , t} and n, m G N, where 

W = max{4 j) (/w) : j G A}. 

For the rest of the proof, we take a point x = G S such that both f |2.5j) and 

fLT6|) hold for x. 

Fix n G N. According to (11. 8p . T~ 1 M Xl ... Xn T is a partitioned matrix of the form 
i<j,k<t, where each J3W> ) is a ci,- x matrix given by 

(2.7) = V Ai j ' yi) A^' y2) ■ ■ ■ A^ n ~ 1,k \ 

i<»l<V2<—<l/n-l<fc 

It is easy to check that the number of words 3/13/2 ■ ■ • 2/ n -i G {1, . . . , t}™" 1 , satisfying 
the restriction j < yi < y2 < ■ ■ • < Vn-i < is bounded above by h(n) = (2n)*. 
Furthermore, each such a word .72/12/2 • • -Vn-ik can be written cLS (X^ 0*2, ' ' ' 

s , where 

s G {1, . . . , £}, j = di < • ■ ■ < a s = k, and ni, . . . , n s G N with n\ + ■ ■ ■ + n s = n + 1. 
Hence 

(2.8) A^ yi) A[ yi ' y2) ■ ■ ■ Aty"- 1 '® = W 1 Ai ai ' a2) W 2 Ai a2 ' a:i) ■ ■ ■ W s ,A^' a ^ W„ 
where 

if rij = 1, 
AB„ 0+ ... +ni _i ir n t > i 

for all % G {1, . . . , s}. Here Idxd is the d x d identity matrix and tlq = 0. Observe that 
(12. 6p gives 

||Wi|| < C{x) exp ((rii - 1) max{W, -A}) exp ((m H h m - l)e) 
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for all i G {1, . . . , s}. Hence, by (12. 8p . we have 

s 

\\a^a^ ,V2) ■ ■ -4r i,fc) n < L *~ l II ll^ll 

i=l 

< L s ~ 1 C{x) s exp((n + 1 - s) max{tV; -N}) exp(nse) 

< DL t C(x) t exp(nmax{W, -N}) exp(nte), 



(2.9) 



where 

L = 1 + max{Pp 2) || : j u j 2 G {1, . . . , t} and i G {1, . . .,£}}, 
D = max{l, exp((i + 1) max{W/, -iV})}. 
Therefore, by (I2.7p - (I2.9I) . we have the estimate 

||T- 1 Af xl ... x „T|| < t 2 max{\\B^\\ : j, k G {1, . . . , t}} 

< t 2 h{n)DL t C{x) t exp(ra(max{W / , -N})) exp(nte) 
for all n G N. Combining this estimate and (12.51) yields 

M,(/i) = lim - log ||T- 1 Af Xl .. <BB r|| < max{W, — iV} + te. 
Letting N — > oo and e — > 0, we get 

M»(/x) < W = max{4^(/i) : j G A}, 
which finishes the proof of part (i) of Theorem 11.71 □ 

3. Extensions and remarks 

For an invertible matrix M G IR dxa! , following [4], we define the singular value 
function of M as 

<P(M) = a x {M) ■ ■ ■ a k {M)a k+l {Mf-\ 

where < q < d, k is the integral part of q, and cXi(M) is the z'-th largest singular 
value of M. For q > d, we put (p q (M) = \ det(M)\ q ^ d . It is known (see [H Lemma 
2.1]) that <f) q is sub-multiplicative in the sense that 

<P q (M 1 M 2 ) < (f> q {M 1 ) ( j)i{M 2 ) 

for any two invertible matrices M l5 M 2 G M dxd . For a given family of invertible 
matrices {Mi} e i=l C IR dx<i , similar to ( 11. lft . we define 

(3.10) P+(q) = lim - log V 9 (Mj). 
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For n G A^o-(E), we define 

(3.11) <FM = lim - V //(J) log^CMj). 

Jes n 

Then by [151 Theorem 2.6], or more generally by (31 Theorem 1.1], we have the 
following variational principle 

P*(q) = max{(f>l(ji) + h(ji) : fi G A4 ff (£)}. 

Similarly we can study the structure of the equilibrium states of P^(q). It is easy 
to see that Theorem 11.71 remains true for P^(q) when 0<q<l or q>d — 1. 
Observe also that it is true when q is an integer: if M Aq is the q-th exterior product 
of M G R dxd (i.e. the ( ) x ( d ) matrix whose entries are the qx q minors of M), then 

ai {M Aq ) = a 1 {M)---a q {M) = 9 (M). 

This gives a partial answer to [HH Question 6.3]. 

Question 3.1. When using ( 13.101) and (13. lip instead of (11. II) and ( 11. 2ft . o?oes some- 
thing like Theorem \l.H\ hold for q G [1, d — 1] \ N? 

We remark that some assumption was given in [5] so that an analogue of (11.51) 
(where || • || is replaced by 9 (-)) holds; and for such case, an analogue of Proposition 
O holds for (cf. Theorem 5.5]). 
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